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Abstract —In this work, we maximize the secrecy rate of the 
wireless-powered untrusted relay network hy jointly designing 
power splitting (PS) ratio and relay heamforming with the 
proposed global optimal algorithm (GOA) and local optimal 
algorithm (LOA). Different from the literature, artificial noise 
(AN) sent by the destination not only degrades the channel 
condition of the eavesdropper to improve the secrecy rate, but 
also becomes a new source of energy powering the untrusted 
relay based on PS. Hence, it is of high economic benefits and 
efficiency to take advantage of AN compared with the literature. 
Simulation results show that LOA can achieve satisfactory 
secrecy rate performance compared with that of GOA, but with 
less computation time. 

I. Introduction 

Recently, simultaneous wireless information and power 
transfer (SWIFT) became an emerging solution for prolonging 
the lifetime of energy-constrained wireless nodes and draw 
significant attention in the cooperative transmission. Since 
SWIFT enables receivers to harvest energy and encode in¬ 
formation from the same wireless signal, it makes the most 
efficient utilization of wireless spectrum for both information 
and energy transfer. SWIFT has drawn a great deal of research 
interests IIl-ll. Two practical schemes, namely power split¬ 
ting (FS) and time switching (TS), were proposed in HI and 
E). For TS, the receivers switch over time between infor¬ 
mation decoding and energy harvesting. FS enables wireless 
nodes to split some power for information decoding and use 
the remaining power for energy harvesting. On the other hand, 
energy harvesting powered relay systems have gained signifi¬ 
cant interests, where the relays with SWIFT use their harvested 
energy to power the forwarding of the sources’ information 
without external power supplies. In JS), time switching based 
relaying (TSR) protocol and power splitting based relaying 
(FSR) protocol were proposed for wireless-powered relays 
systems. The authors in m investigated energy harvesting 
relays networks where relays switched over time between 
harvesting energy and volunteering to forward the transmitter’s 
information. In Q, a multiple-antenna relay system with 
SWIFT was studied and a “harvest-and-forward” strategy was 
proposed to maximize the achievable rate. The authors of il 
investigated the beamforming design with simultaneous energy 
harvesting to improve the max-min signal-to-interference-plus- 
noise ratio (SINR) in a multiuser two-way multi-antenna relay 


systems. 

On the other hand, with the development of wireless net¬ 
works, the issues of privacy and security have attracted much 
interest and attention due to the broadcast nature of the wire¬ 
less medium. The cryptographic approach using secret keys 
at upper layers has been traditionally used for the security of 
communication systems. Due to the cost of the key distribution 
and key management, information-theoretic security at the 
physical layer became an important complementary to the 
traditional cryptographic approach. Thus, substantial research 
efforts have been dedicated to information-theoretic physical 
layer security a-nsi, the concept of which was defined 
by Wyner cni. By exploiting spatial diversity, a relay may 
also provide the secrecy capacity by assisting the source- 
destination transmission or acting as a jammer (e.g. CD). 
However, another security threat arises since the information 
transmitted by the source may be eavesdropped by an untrusted 
relay ifT^ . Interestingly, in lfT2l . it was shown that the secrecy 
rate may be enhanced even with the help of an untrusted relay. 
In m, the authors studied joint secure beamforming for an 
AF untrusted relay system. The MIMO AF untrusted relay 
beamforming was further extended to a two-way untrusted 
relay system in US. To further improve the security of the 
source-destination transmission, artificial noise (AN) (sent by 
jammer ifTTI . source iflTl or destination ifTSll ) was proposed in 
m as another important solution for physical layer security. 
Traditionally, AN is only used to degrade the channel con¬ 
ditions of eavesdroppers, hence, it may not be efficient and 
economic from the aspect of green communication. 

More recently, a handful of works studied the SWIFT with 
the consideration of physical layer security QD-Illll. The 
authors in ll2Tll studied secure relay beamforming for SWIFT 
in one-way relay systems with an external energy-harvested 
receiver. In ll22l . the authors considered a wireless-powered 
jammer system, where the jammer can harvest energy from the 
wireless signal and use it to interfere with the eavesdropper. 
In the present literature, the energy of AF relay, especially 
the untrusted relay, is considered mostly to be supplied by the 
power grid or battery. To the best of our knowledge, secrecy 
issues in the wireless-powered untrusted relay systems have 
not been studied yet, which motivates this work. And it seems 
reasonable and intuitive to realize the energy harvest at the 


untxusted AF relay node to make full use of AN. 

In this work, we consider a wireless-power untrusted relay 
network and maximize its secrecy rate by jointly designing 
PS ratio and relay beamforming with the proposed global 
optimal algorithm (GOA) and local optimal algorithm (LOA). 
The main difference from the literature is that AN sent by 
the destination not only degrades the channel condition of the 
eavesdropper to improve the secrecy rate, but also becomes a 
new source of energy powering the PS based untrusted relay. 
Hence, it is of high economic benehts and efficiency to take 
advantage of AN compared with the literature. 

Notation: We adopt the notation of using boldface for 
column vectors (lower case), and matrices (upper case). The 
hermitian transpose is denoted by the symbol j-)^. For a 
complex scalar x, its complex conjugate is denoted by x*. 
E[-] and CAf{-) denote the statistical expectation and complex 
Gaussian distributions, respectively. For matrix X, vec(X), 
X~^, X'^, and Tr(X) represent the vectorization, inversion, 
transpose and trace of matrix X, respectively, and X ^ Y 
stands for the kronecker product of X and Y. For a vector 
X, we use ||a ;||2 to indicate its £2 norm. And I is the identity 
matrix with corresponding dimensions. 

II. System model 

Consider a scenario with three nodes: an untrusted AF relay 
(R), source S and destination D, where only R is equipped 
with multiple antennas, and the scenario with single-antenna 
R is regarded as a special case in this paper. Due to the space 
limitation, the multi-antenna case where all the nodes own 
multiple antennas is considered in the journal version. 

In this paper, the considered relay R is untrusted and 
assumed to eavesdrop the conhdential signal of D when it 
helps the SD transmission, and it can be treated as a legitimate 
user ifT^ - lfTSl with different service from that of D. In FiglT] 
at the first hop, S transmits the confidential signal to R, 
meanwhile, D sends AN to R to degrade the eavesdropping 
channel of R. In contrast to the literature, R can not only 
process the received signal, but also harvest energy from it as 
the source of power for the next-hop transmission at the same 
time. At the second hop, R transmits the processed signal to 
D depending on the harvested energy. Notice that both R and 
D are half-duplex, i.e., R and D can not receive and send 
signal, simultaneously. Therefore, SD link is not under our 
consideration. 

Denote the channel vectors from S to R, from D to R, and 
from R to D by hsr S h^r S and hj-d S 

(j^AfrXi, wjjere Nr is the number of antennas at R. At the 
first hop, S transmits the signal to R, meanwhile D sends AN 
to R with the purpose of confounding the eavesdropper. The 
received signal at R is given by 

Vr — ^sr^s ^dr^d Llj., ( 1 ) 

where G represents the additive white Gaussian noise 
(AWGN) associated with relay node following the distribution 
CA£(0,cr^J), Xs represents the secure information sent by S 
to D, meanwhile, Xd is AN sent by D to degrade the channel 
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condition of R. And Xs S C and Xd G C represent the transmit 
symbol of S and D with powers E[xsxl] = Pg and Elxdx’^] = 
Pd- p denotes the ratio of power split for the energy harvesting, 
and the information processing is (1 — p) at the relay node. 
Hence, the parts of received signal for energy harvesting and 
information processing are written, respectively, as follows 

Vr — v/P i^sr^s ^^dr^d Ltr) ; (2) 

Vr — y/l P i^hgrXs hdrXd Xlr') -f TIq. (3) 

where tie G C^'’ is the additive noise vector introduced by 
signal conversion from the passband to baseband with the 
distribution CAf{0, cr'^I). 

Based on ||23], we denote A{p) := (1 — p)Pdhdrh\^ + (1 ~ 
p)o-rI -f (Jr I- Furthermore, we can rewrite Q as 

-tJr = \/l - pA~^{p)hgrXs + 

A~^{p) phdrXd + v/1 - prir + ,(4) 

and the information rate achieved at the untrusted relay can 
be given by 

Rr{p) = log 2 (^1 -f (1 - p)PshlrA~^{p)hsr'^ - (5) 

Resorting to the PS scheme, we can briefly write the 
expression of harvested energy at R as 

-P™ = V [pPs \\hsr\\l + pPd \\hdr\\l + per^Tr (/)|, (6) 

where rj G (0,1] denotes the energy conversion efficiency from 
signal power to circuit power. 

Different from simply applying the amplified scalar a to 
tackle the received signal at the single-antenna R in the 
Remark [H in this paper, we adopt F as the processing matrix 
to deal with the received signal and forward it to D at relay 
node. Subtracting the AN term — phl^FhdrXd sent by 
itself, the received signal at D can be given by 

yd = hl^F (^ a /1 - phsrXs + ptir + + Ud, ( 7 ) 

where Ud is the AGWN introduced by the receiver antenna at 
D with the distribution CAf{0,a^). For simplicity of descrip¬ 
tion, we make cr^ = = tr^ in the following parts. 
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III. Problem Formulation For Secrecy Rate 
Maximization 


To begin with our problem formulation, we will provide a 
brief background on the introduction of Secrecy Rate im- 
M- In the presence of the eavesdropper, the secrecy rate is 
the information rate that the source can securely transmit to 
the legitimate user. In this section, we denote Rsr{F,p) as 
the secrecy rate, given by 

R^riF, p) = i [i?d(F, p) - RriF, p)]+ , (9) 

where [a]"*' = max(0,a), and i is from the fact that the 
source-destination transmission takes place in two time slots. 

For the MIMO relay case, the optimization problem can be 
written as 


max 

P,F 

s.t. 


Rsr{F,p) 

il-p)Ps \\Fhsr\\l + il - p)Pd\\Fha^\l + 

{l-p)a^,\\F\\l+al\\F\\l<Pr, 


( 10 ) 


which can be briefly rewritten as 


max 

pj 

s.t. 


A+__V_^ 

V /tC(p)/ + pEf + Vl + (1 - p)PshtA-Hp)h,J 

/lT(p)/<P™^ (11) 


where 


/ 

B{p) 

Cip) 

E 

Gip) 

Jip) 

Tip) 


vec(F), 

Psil-p) 

(l-p)[( 

:( 




a-1® hrdhl^ 


^Id 


(1 - p)Ps 

(hsrhl^'j 0 I 

(1 - p)pd 

(hdrhl^'j 0 I 


{Gip) + J{p) + il-p)a^If + a^Jf), 


where If € has the same dimension as G{p), and 

the properties, such as vec{AXB) = {B'^ 0 A)vec(X) and 
Tr(X^X 2 ) = vec(Xi)^vec(X 2 ), are applied to the above 
equalities. Note that the energy constraint in ( fTOl i means that 
the considered relay is wireless-powered, i.e., the relay is 
powered by its harvested energy without any external energy 
supplies. 


Remark 1: When S, D and R are equipped with the single 
antenna, the achievable rate of R and D are given by 


Rr(p) = log 2 1 -f 


(l-p)P,lh, 


Rd(p) =log2 1 -f 


(1 - p)Pdlhdrj^ + (1 - P)^d + 

a^il - p)Ps\hljlsr\^ 


a2(l - p)\hrd\^(rr + a^\hrd\^al + aj 


t2 1 ’ 


where a is the amplified factor with the consideration of power 

I p™ 

constraint at the relay, a = ^ • 

And the secrecy rate is 

R,rip)=\[Rdip)-Rrip)\^ . ( 12 ) 

which is a function of p in the range of (0,1), and can be 
perfect solved by setting = 0 (resorting to Matlab), 

where is the derivative of Rerip) '''l^h respect to p. Note 
that there might be more than one point making the derivative 
equals zero, and we treat the one yielding the maximal 
secrecy rate Rgrip*) as the optimal PS ratio p*. If there is no 
point at which the derivative is zero, the secrecy rate is zero. 


A. Global Optimal Solution 

It is obvious that (fTTT i is non-convex and difficult to solve 
with the two variables (p, /), simultaneously. In this sub¬ 
section, we design the global optimal algorithm (GOA) for 
problem (fTTT i with one dimension search of p in the interval 
of [0,1]. In each round, the following optimization problem 
with the given p is considered to obtain the optimal solution 

r, 


max 

/ 

s.t. 


f^Bip)f 

fCip)f + fEf + a^ 
f^Tip)f < P™. 


(13) 


However, the above optimization problem is still a non-convex 
problem. Hence, we resort to semidehnite relaxation (SDR) to 
drop the rank-one constraint (rank(JA) = rank(//l) = 1), and 
obtain the relaxed problem as follows 


TriBip)X) 

x"" Tr(C(p)X) -f TriEX) + 

s.t. Tr(T(p)X) < 

X P 0, 


(14) 


which is a quasi-convex problem. Thus, we adopt the Charnes- 

Cooper transformation to make the quasi-convex problem (fT4l) 

7 

into a convex SDR To this end, we set X = f, e > 0, and 
recast (Ell as 

max Tr(B(p)Z) (15) 

Z,(, 

s.t. Tr(C(p)Z) -f TxiEZ) -f = 1, 
Tr(T(p)Z) < ^P™, 

Z P 0,C > 0, 
































the optimal solution of which can be efficiently obtained by 
solvers, e.g., CVX. Denote as the optimal solution 

of problem (fTSl l with given p, and we achieve the optimal 

F* 


solution of X* by = -jf. Note that the optimal X with 


rank-one property, is the sufficient and necessary condition for 
the equivalence of problem (fOl l and (fT4l) . which is guaranteed 
by the following lemma: 

Lemma 3.1; ll24l Theorem 2.1] Suppose that T € is a 
complex Hermitian positive semidefinite matrix of rank r, and 
Ai, A 2 € TL^ be two given Hermitian positive semidefinite 
matrices. Then, there is a rank-one decomposition of T. 


i=i 


such that 

■v]AiVj=^ 


Tr(AiT) 


A 2 VJ = 


Tr(A2T) 


(16) 


= ,r. (17) 


In our problem, if we want to guarantee the rank-one 
property, the following equalities should be satisfied; 

TT{C{p)Z)+TiEZ)+^a^=TT(c{p)ff^)+TT^ff^+^a^m 


Tr{T{p)X) = Tr{Rip)ff 


(19) 


Based on Lemm^3T] we can recover the rank-one solution 
f* from X*p = f*pfy' with given p. For the sake of brevity, 
denote /(p,/) = /i(p,/)/ 2 [f(p)], where hip, f) = 1 + 

l+{^—p)Pshlj.A~^ ip)hsr- And the optimal value of objective 
function with given p is denoted as f{p,fp). Hence, the 
optimal solution (p*, /*) to problem (fTTT i can be obtained by 


(P*, /*) = arg max /(p, fp). (20) 

pe[o.i] ^ 

The details of GOA is summarized in Algorithm 1 as 

Algorithm 1: The global optimal algorithm 
0: Define Ap = = e. where M is a large constant. 

1: For ) = 1, • • • , M, do Step S0^S2, 

and seek for the global optimal solution (p*, /*) for (fTTI) . 
SO: Denote p = jAp. Solve problem (fTSi) with the given 
p, and denote Zp as the optimal solution to (fTsl i: 

SI: Obtain X* by X* = Do the following 
procedures to obtain the rank-one solution /* 

If rank(Xp=l 

The optimal value of problem (fOT l is /p 
Else 

Resort to Lemma (13.11) to seek for a rank-one 
decomposition of Xp i.e., X7 = zz'^ and /7 = 2 :, 
which satisfies the equalities of ( fTSl l and ( fT^ . 
f*p is the optimal solution with given p. 

End If 

S2: Remark and /(p, /p as the optimal solution 
and objective value for given p, respectively. 

2: Choose the optimal solution from the following equation 

iP*,f*) = argmaXpg[o,i] fjpjp). 


where e is the given tolerant error. 


B. Local Optimal Solution 

Although the global optimal solution (/*, p*) can be found 
as described in Section IIII-AI the high computation time of 
one dimension search in p is expensive. Hence, we want to 
seek for some other algorithms which can better balance the 
performance of system and the cost of computation. In this 
section, a well-designed two-variable block coordinate decent 
(BCD) method is introduced to solve problem (fTTT i converging 
to a local optimal solution, which is proved in Il25l Theorem 
4.1]. 

As it is mentioned above, cr^ = tr^ = = cr^, then 


A-\p) 


(1 - p)Pdhdrh\^ + (2 - p)cr^I 


In order to get more insight of p, we make some transformation 
of term A~^(p) to present it as the function of p instead of 
matrix formulation, like 


.-l, {2-p)a^^<i^dr^drI Pdhdrh\,, 

(2 - p)a^ + (1 - p)Pdhlhdr ’ ^ ^ 

the details of this procedure will be shown in Appendix 
It is noted that different from the global optimal algorithm, 
searching p* to yield the maximal secrecy rate in the interval 
[0,1], the local optimal search for p implies that we should 
select p making the secrecy rate greater than zero to guarantee 
the secure source-destination transmission. Notice that the 
choice of p decides the nature of secure transmission in this 
relay networks. For example, small p means that there is more 
information eavesdropped by the untrusted relay (see 0), 
and less harvested energy for the second-hop transmission, 
resulting in lower achievable rate for D. Meanwhile, although 
the large p (approaches to 1) leads to more harvested energy 
and less information leakage, the confidential information is 
weakened and difficult to be distinguished from the noise at 
the destination. 

Therefore, it is important to choose a suitable p for this relay 
networks. In this subsection, we resort to BCD to update p and 
/ in sequence. Due to the property of BCD, the initial point 
has a great effect on the solution of problem (ITTI) . Hence, we 
choose multiple initial points, such as pp and , j = 1, • • • , J, 
to achieve different objective values, and the maximal one is 
set to be the local optimal value. 

For each initial point (p, f) (denote p = p’ and f = f\ Vj 
for simplicity), problem (ITTIi with a given fixed to achieve 
the temporary optimal solution p\ in the fc-th iteration, can be 
written as 


max 

p 

s.t. 


fi+__V_^ 

fiTip)fp < P™. ( 22 ) 


Plugging (I 2 TI 1 in (l22l i. it is easy to verify that the objective 
function of problem (l22T i can be rewritten as an univariate 
function of scaler p with high order term, and the constraint 
with respect to p is convex. Therefore, the optimal solution 
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p*f. of the above problem can be obtained. We check the 
monotonicity of with respect to p (the details are 

given in Appendix IbI). choose the points where the derivatives 
of /(p, fj.) with respect to p equal to zero, calculate the values 
of the objective function in (l22l) with those chosen points, and 
select the one with the maximal value. If there is no point 
making the derivative of f{p,fk) with respect to p equal to 
zero in the range of [0,1], it means that /(p, f^.) is monotonic, 
and the point at the boundary of the constraint will be chosen 
to be the optimal one. The reason is given in the Remark |2 

Remark 2: As / is a smooth function, V/ should be a 
continuous function. We have that if Vp/ can not obtain 0, 
Vp/ must be positive or negative in the whole interval [0,1]. 
If not, let us assume that Vp/(pi, ff.) and Vp/(p 2 , f^) have 
different signs. Because of the continuity of Vp/, there should 
be a point pa e (pi, P 2 ) such that Vp/(p 3 , fk) = 0. This is a 
contradiction. 

Then we can update / depending on the obtained pfc+i 
(where pfe+i := p\) in the (fc+ l)-th iteration through solving 
the following relaxed optimization problem referred to the 
procedures of ([T4l) to dfSl l 

max Tr(JB(pfe+i)Zfc+i) (23) 

s.t. Tr(C(pfc+i)Zfc+i) + Tx{EZk+i) + = 1, 

Tr(T(pfe+i)Zfc+i) < 

the optimal solution of which can be efficiently obtained by 
solvers, e.g., CVX. Denote as the optimal solution 

of problem (|2^ in the [k + l)-th iteration, and we achieve 

Z* 

the optimal solution of X* by . Resorting 

to lemma 13.11 we can get a rank-one decomposition from 
= n+Jllx, i.e., n+i is the optimal solution in the 
[k + l)-th iteration. 

In the following, we present the procedures of the local 
optimal algorithm (LOA) in Algorithm 2. 


Algorithm 2: The LOA based on BCD method 
0 (Initialization); Randomly generate feasible points(p-^ , 

V/, and let the initial objective values be f{p^, 

Do Stepl-2 for each j, denote p = p^ and f = 
for simplicity, and set A: = 0. 

1 (Block Coordinate Maximization): 

(1) Solve Problem (l22l i with f{p, /j,), and 
get optimal solution pfc+i; 

(2) Solve Problem (l23t with f{pk+i, /), and 

get optimal solution resorting to Lemma (TTl 

Set A = |/(pfe+i,/fc+i) - f {Pk,fk)\/\fiPk+i,fk+i)\- 

2 (Stopping Criterion); If max(0, A) < e, stop and return 
(pfe, ff.)-, otherwise, set k := k + 1 and go to Step 1. 

3 (Local Optimal Solution): Denote {pi, fl.) and f{pl.,fl,) 
as the solution and objective value of problem (fTTI i with 
the given initial point (p^ , /^), hence the local optimal 
solution is chosen by 

(p*, /*) = argmaxj /{pj, fj). 



Fig. 2. Verification of Remark [T] that the optimal p* is obtained when the 
derivative (p*) equals zero, Ps = Pd = 40dBm. 

IV. Simulation Results 

In this section, we investigate the performance of the 
proposed algorithms via numerical simulations. The path loss 
model for the energy harvesting relay channel is Rayleigh 
distributed and denoted by where |/3| and d represent 

the short-term channel fading and the distance between two 
nodes (S and R, D and R, R and D). |/3p follows the 
exponential distribution with unit mean. We set the noise 
powers as cr^ = cr^ = cr^ = cr^ = OdBm, Vi, the number of 
antennas at R as 2, and the energy conversion efficiency p = 1. 

To get more insight of Remark [T] we verify the conclusion 
of Remark [T] using one channel realization. For the single¬ 
antenna untrusted relay case, the optimal PS ratio p* can be 
achieved by setting = 0. As shown in Fig|2l the point 
with the maximal secrecy rate in the upper subfigure corre¬ 
sponds to the point p* where = 0 on the horizontal 

axis in the second subfigure. In addition, we conclude that both 
small p and large p (approaches to 1) yield smaller secrecy rate. 

The secrecy rate performances for GOA and LOA in MIMO 
untrusted relay networks are given in Fig|3] In the low SNR 
region, the performance achieved by LOA quite approximates 
the one in GOA. The reason is that the key element affecting 
the performance in that SNR region is the harvested energy 
powering the next-hop transmission. Although we can get the 
optimal p* resorting to GOA, it just helps a little. 

In contrast, the dominant factor for the performance is the 
power splitting ratio p in the high SNR region. It is the fact that 
the higher the SNR, the larger gap between GOA and LOA for 
the fixed AN power in Fig 13 Interestingly, as the power of AN 
{Pd) increases, the performance of LOA is getting much closer 
to that of GOA for the whole region of SNR. For instance, at 
SNR=50dBm, the relative performance ratio is growing from 
79.6% to 86.3% as Pd raises from 40dBm to 50dBm. 

The comparison of computation time for GOA and LOA 
is shown in TableH] with different given e. The results are 
averaged among 100 channel realizations, and the unit of 
measurement for computation time is secoiid{s). It is obvious 
that LOA takes much less computation time than GOA. We 





















Fig. 3. Comparison of the secrecy rate achieved by different algorithms. 



Fig. 4. Comparison of relative power used (relative to the power consumption 
of EPR case with = 40dBm and Pr = 43dBm. 


power Pr = 43dBm is almost twice times than those of 
GOA and LOA in the low region of SNR. Because WPR 
can not harvest enough energy to achieve better performance 
in that low SNR region compared with EPR. By contrast, in 
the high SNR region, WPR can harvest sufficient energy to 
support its next-hop transmission and obtain higher secrecy 
rate performance with reasonable PS ratio for the information 
processing, but consumes less energy compared with EPR 
in EiglD As for EPR, more information is eavesdropped by 
the untrusted relay as the source power Pg is increasing in 
the high SNR region, hence, its secrecy rate performance is 
reducing. Prom the aspect of secrecy rate performance, WPR 
has advantage over EPR for the reason that it can harvest 
energy from the wireless signal and adjust the PS ratio to 
achieve better performance with less information leakage and 
power consumption, especially in the high SNR region. 

V. Conclusion 

In this paper, we study the secrecy rate maximization 
problem in the wireless-powered untrusted relay networks 
according to joint power splitting and secure beamforming 
design. We propose two algorithm named as Global optimal 
algorithm (GOA) and Local optimal algorithm (LOA) with 
well established convergence result to jointly optimize the PS 
ratio and relay beamforming, and LOA can achieve satisfac¬ 
tory secrecy rate performance compared to GOA, but with 
less computation time. We also show that the wireless-powered 
untrusted relay has advantage over external-powered untrusted 
relay at the aspect of secrecy rate performance but with less 
power consumption, especially in the high SNR region. 


TABLE I 

Comparison of The Computation Time in GOA and LOA for per 

CHANNEL REALIZATION. 


Appendix A 

Transformation of A~^{p) 


Method 

e 

10"-“ 

10"'" 

10-“ 

Local Optimal Algorithm (LOA) 

0.2896 

2.9868 

7.4S35 

Global Optimal Algorithm (GOA) 

0.9706 

6.1656 

47.9588 


We rewrite A ^{p) to represent it as the function of p 
instead of matrix formulation, as follows 


A-\p) 


(2 - p)cr^I + (1 - p)Pdhdrh\^ 


can figure out that LOA is relatively efficient, and it can 
achieve satisfactory secrecy rate performance compared to 
GOA, especially when Pd is large, but with less computation 
time. 

To find more knowledge of the benefits of the wireless- 
powered relay, we compare the performances between the 
wireless-powered relay (WPR) and the external-powered relay 
(EPR, i.e., the power of the relay is supplied by the power grid 
or battery, rather than harvested from its received signal) in 
Pigl3] We also show the relative power consumption of the 
three algorithm (GOA, LOA and EPR) in EiglD and denote 
GOA and LOA as WPR-GOA and WPR-LOA, respectively. 
It is clear that WPR consumes less power than EPR because 
the power of WPR is supplied by energy harvesting from the 
received signal, rather than the external power for EPR case. 

The secrecy rate of EPR with the given constant relay 
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where the second equation comes form the definition of 
Sherman-Morrison-Woodbury formula ll^ . which is 
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Appendix B 

The Derivative of f{p, f) with respect to p 

Recall the those denotations in Section IIII-AI of /(p, /), 
fi{p, f) and f 2 \t{p)], we have the partial derivative of fi with 



























































respect to p like 


-P f 

dhipJ) 


\h,rhl)®{hrdhl^ f(^f^Ef+a^ 


(fc{p)f + fEf+a^' 

And the derivative of /2 with respect to p, is given 
dMt{p)] _ 1 dt{p) 
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where t{p), taking account of the new transformation of 
A~^{p), is written as 

t{p) = 1 + mn, 


where m = 




{2-p)a^ 


i^-pW 


and n = 


(i-p)A 


{2-p)o^ + {l-p)Pdhl^hdr- ' 

The derivative of f{p, /j.) with respect to p is given by 
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